A geometric approach to a generalized virial theorem 



a* 

00 



Jose F. Carinena a \ Fernando Falceto^ and Manuel F. Rariada c ) 

Departamento de Fisica Teorica, Facultad de Ciencias 
Universidad de Zaragoza, 50009 Zaragoza, Spain 

September 21, 2012 



i-^ ■ Abstract 
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Lagrangian systems. The possibility of establishing virial's like theorems from one- 
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1 Introduction 



In 1870 (see [Tj for an historical account), Clausius introduced the virial function for a one 
particle system 

G(r, v) = m r ■ v (1) 

for studying the motion of a particle of mass m under the action of a force F. Then, using 
such a function and the Newton second law, he proved that when either the motion is 
periodic of period T, or the motion is not periodic but the possible values of the function 
G are bounded and we take the limit of T going to infinity, the time average over a time 
interval T of kinetic energy E c (v) is given by 

((£ c (v)» = -i«r-F», (2) 

so that, in the particular case of a conservative force, he obtained 

»)))^({r-VF)>, (3) 

If the potential V is homogeneous of degree k, Euler's theorem of homogeneous functions 
implies that r- W = k V, and therefore, 2 ((E c (v))) = k ((V(r))), that leads to the following 
values for the averages of the kinetic energy and the potential 

«*(v)» = j£, «VW» = j£, 

where E is the total energy [2]. For instance in the harmonic oscillator case, with k = 2, 
and in the Kepler problem with negative energy, with k — — 1, we obtain, respectively, 

((E e (v))) = «y{r))) = ±E, 
«£ c (v))) = -E, ((V(v))) = 2E. 

The virial theorem was originally introduced in classical statistical mechanics (that was 
the original matter studied by Clausius) but later became important in many other differ- 
ent branches of physics. The quantum mechanical version of the theorem is due to Born, 
Heisenberg, and Jordan [3] and presently it is a tool frequently used in many body quan- 
tum mechanics (e.g. systems of fermions) and in molecular physics; in addition, this matter 
has also been related to certain fundamental questions arising in nonrelativistic quantum 
mechanics as the Hellmann-Feynman theorem (see [1]-[IH] and references therein). We also 
mention the importance of the virial theorem in solar and stellar astrophysics [15] . The 
important point therefore is the wide range of applicability of the virial theorem with ap- 
plications ranging from dynamical (even relativistic) and thermodynamical systems, to the 
dust and gas of interstellar space, as well as cosmological considerations of the universe as a 
whole and in other discussions concerning the stability of clusters, galaxies and clusters of 



2 



galaxies [20J , or giving information on the masses of bound systems, being the main reason 
we think that dark matter exists J2TJ [22] . Of course this virial theorem provides less infor- 
mation that the equations of motion but it is simpler to apply and give us some information 
on systems whose complete analysis may defy description. 

After this standard approach to virial theorem several questions appear. For instance, the 
role played by the function G and also the reason why the relation is simpler for power law 
potentials. There are more general results known as hypervirial theorems introduced in [5] 
whose natural framework is the theory of locally Hamiltonian dynamical systems. Moreover, 
quantum mechanics is also a particular case of this more general theory which suggests 
us to analyse the problem from this new perspective, that then enables us to develop the 
corresponding theorem in quantum mechanics. The main point is that the standard theory 
is always related with scaling and dilation transformations [TQ1 EED S [231 [21]; there is no 
room for such transformations when the configuration space is no longer a linear space but 
there exist, for example, holonomic constraints. 

Concerning the last question, the search for generalizations is a matter that has been 
studied by different authors from different viewpoints (in fact the quantum mechanical virial 
for the expectation values of a quantized system can be considered as a generalization of 
the original virial theorem). Here we mention a version for discrete maps [22] or a virial 
theorem in a spherical geometry [25]. An interesting result is that it has been related with 
the Noether's theorem and with the invariance properties of the Lagrange function [271 121] ■ 

The main aim of this paper is to develop a deeper analysis of the virial theorem, both 
in classical and in quantum mechanics, using the modern theory of Hamiltonian systems in 
symplectic manifolds as an approach, and therefore including Lagrangian formulations. 

The paper is organised as follows. In Section 2, a generalized virial theorem is established 
in the framework of the theory of Hamiltonian systems in symplectic manifolds. In Section 
3, the virial theorem is related with the Lagrangian formalism by considering the Lagrangian 
systems as Hamiltonian systems in tangent bundle manifolds; we prove that not only point 
symmetries of L are useful in classical mechanics but more general transformations modifying 
the Lagrangian by a factor, can be used to establish virial type relations. Two generalizations 
of the virial theorem are studied in Sections 4 and 5. In Section 4 we consider the case of 
systems with a position dependent mass and in Section 5 we relate the theorem with the 
theory of non-strictly canonical transformations. We study in Section 6 the virial theorem 
in quantum mechanics using quantum symplectic formalism as an approach. Finally, in Sec. 
7 we make some comments and present some open questions. 

2 Virial theorem for Hamiltonian systems in symplec- 
tic manifolds 

To make the paper self-contained and to introduce the notation, we recall some concepts of 
the geometrical approach to the theory of Hamiltonian systems. Then we consider the virial 
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theorem from a geometric perspective. 

If (M, u) is a symplectic manifold, the Hamiltonian vector field Xp defined by the function 
F on M is the solution of the equation 

i(X F ) u = dF 

A dynamics is given by choosing a Hamiltonian H through the corresponding Hamiltonian 
vector field Xh- 

For instance, the phase space M of the Hamiltonian systems studied in classical mechanics 
is the cotangent bundle T*Q of the configuration space Q [29] -[32]. They are endowed with 
a canonical exact symplectic structure uo that in cotangent bundle coordinates (qi,Pi), i = 
1, . . . , n, is given by uo = dqiAdpi, i.e. ujq = —d9 , with 9q = Pi dqt, where the summation on 
repeated indices is understood. The Lagrangian formulation for regular Lagrangian functions 
L is also an example. The manifold is the tangent bundle M = TQ and the symplectic 
structure is not canonical but depends on the choice of the function L. 

The Poisson bracket of two functions F\ and F 2 in a symplectic manifold (M, u) is defined 
as the symplectic product of the corresponding Hamiltonian vector fields 

{Fi, F 2 } = u>o(Xp 1 , Xp 2 ) = Xp 2 Fi . 

It can be verified that, when written in Darboux coordinates, the integral curves of Xjj 
are solutions of Hamilton equation and for the Poisson brackets we recover the expression 
usually given in classical mechanics. 

The flow of the Hamiltonian vector field Xjj on a symplectic manifold (M, u) commutes 
with the action of X^] thus, if F is a function in M we have 

4>*{X H F) = X H {<jf t F) = 1[0*(0*F)]| S=O = ±(<p* s+t F) ls=0 = ^[KF] lu=t . 

For a given function G which is going to play the role of the function (PQ) we obtain 

j t (<P* t G) = <p;({G,H}) = -<p* t (X G H) , 

and integrating this relation with respect to the time from t = to t = T, we arrive at 

i[G o T - G] = -i j\x G H) o4> t dt=± £{G, H} o 4> t dt. (4) 

This result is sometimes known as the hypervirial theorem [S]: If either the motion is 
periodic of period T, or when G remains bounded in its time evolution, taking the limit 
T — > oo: 

= 0. (5) 
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Actually a family of hypervirial theorems was introduced in [5] with an especial emphasis 
on homogeneous first-degree functions in momenta, but the formulation is still valid for any 
choice of the function G. 

In the particular case Q = M 3 the phase space is the cotangent bundle M = T*M. 3 which is 
endowed with its canonical symplectic structure uq. Then the Hamiltonian vector field Xq 
of the Clausius function G given by G(r, p) = r • p, takes the form 




and represents the generator of dilations in M = T*R 3 : the infinitesimal dilations in Q = M. 3 

3 

are generated by D = d/dx 1 and its cotangent lift to the phase space as an infinitesimal 

i=i 

generator of point transformations, gives us Xq. 
When if is a natural Hamiltonian (kinetic term plus a potential) 

H{r, p) = ^ p 2 + V(r) = iZo(p) + V{r), 

the action of Xq on H becomes (XgH)(t, p) = — 2Ho(p) + r • VV. Therefore, the relation 
© is 

i[Go0 T -G] = (<2#o>>-<<r- W>>, 

and, in a periodic motion of period T, or in the limit T — > oo when G remains bounded in 
its time evolution, fl5]) reduces to the standard result in the Hamiltonian framework 

«2ii » = <<r- W». 

Summarizing, the standard virial and hypervirial theorems can be considered as particular 
cases of more general properties stated for Hamiltonian systems in symplectic manifolds. 



3 The virial theorem in the Lagrangian formalism 

It has been shown (see [33l [31]) that a (regular) Lagrangian system is a particular case of 
a Hamiltonian system on a symplectic manifold, on the tangent bundle TQ of the configu- 
ration space Q and with a symplectic structure depending on the Lagrangian function. We 
can therefore translate the theory developed in the Hamltonian formalism to a Lagrangian 
framework. 

Two important geometric ingredients are the Liouville vector field A, that is the generator 
of dilations along the fibres, and the vertical endomorphism S. Given a differentiable function 
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L on TQ, we can construct a semibasic 1-form 9 L e f\ (TQ), an exact two- form ojl E 
/\ 2 (TQ) and an Energy function by 

9 L = S*(dL), u L = -d6 L , E L — A(L) — L . 

If the Lagrangian L is regular then ojl is symplectic and the Lagrangian dynamics is given 
by the uniquely determined vector field X L solution of the equation 

i(X l )uj l = dE L . 

X L satisfies the second-order property S(X L ) = A, and the curves on Q that are a projection 
of the integrals curves of Xl in TQ satisfy the second-order Euler-Lagrange equations. 

To summarize, a regular Lagrangian L determines a symplectic structure ool in TQ and 
the Lagrangian formalism is an instance of the theory of Hamiltonian dynamical systems. 
The virial theorem reduces in this case to 

I[G o T - G] = -I j\x G {E L )) o<p t dt=^ £{G, E L } L o <p t dt, (6) 

where Xq is the vector field such that i(Xc)^L = dG, and {G, El}l = ujl(Xg, Xl). If the 
motion is periodic of period T, or if the function G remains bounded in its time evolution, 
when taking the limit of T going to infinity we obtain 

((X G (E L ))) = 0. (7) 

Let us consider the simple case of a natural Lagrangian defined in Q = M. 3 

L(r,v) = |mv 2 - V(r) . 

The energy function and the symplectic form are given by 

E L (r, v) = | m v 2 + V(r), u L = mdrAdv, 

and when G is the observable function G(r, v) = mr ■ v, then 

*° = £(* , 5F-«'5?). 

i=i x ' 

or written in a simpler way 

X G (r,v) = r- V r -v- V v . 

Note that here the minus sign shows that X G is the Hamiltonian vector field defined by G 
but it is not the lift to the tangent bundle TQ of the infinitesimal generator of dilations in 
Q = IR 3 , but X G depends on the function L. Then we have 

(X G (E L ))(r,v) = r • W(r) - v • (m v) , 
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and we thus recover the original virial theorem ([3]). 

One-parameter groups of point symmetries of the Lagrangian lead to constants of motion 
but infinitesimal transformations that are not symmetries of L can also play a role in es- 
tablishing virial theorems. Next we give a geometric approach to a result of [27], whivh is 
a particular case of a more general result to be given in next section where one-parameter 
groups of non-strictly canonical transformations are to be considered. 

It is important to remember that that if is a diffeomorphism in TQ that is obtained 
from a diffeomorphism ip in the base Q, that is <fi = Ttp, then the following two properties 
[31] are satisfied 

cf)*0 L = Q^l , 4>*E L = E^ L . 

Correspondingly, at the infinitesimal level, for a vector field X G 3C(TQ) that is a complete 
lift, X = Y c , of a vector field in the base, Y G X(Q), we have 

£-x9l = ®x(l) > X{E L ) = E X (l), 

where Cx denotes de Lie derivative with respect to X. 

In the following, given a 1-form a G /^(Q) we denote by a a function on TQ that is linear 
in the fibres defined by a(q,v) = (a q ,v). Then, as 9s is the pull-back of a and Es = 0, 
V = L + a defines the same Hamiltonian system on TQ as L if (and only if) a is closed. 
Hence a vector field X = Y c , Y G £(Q), such that X(L) = a, with a a closed form in Q, 
is a symmetry of the Hamiltonian dynamical system (TQ,u>l, El) defined by L. We can at 
least locally write a = dh where h is the pullback through the tangent bundle projection 
of a function h in Q. Note that, making use of this notation, the property dh = T(h), is 
true for any vector field F on TQ satisfying the second order differential equation condition, 
S(T) = A. 

Theorem 1 Let L be a regular Lagrangian and X = Y c a vector field on TQ such that 
X(L) = aL + dh, with aGl. Then: 

(i) the vector field X is a symmetry of the dynamical vector field Xl; 

(ii) the function G = i{X)9l — h is such that Xl(G) = aL. 

Proof, (i) First, we have 9x(l) — a9 L + dh, and therefore ujx(l) = clojl. Furthermore, 
the energy Ex{l) is given by Ex{l) = a El. Consequently C-x^l = wx(i) = clcol and 
CxEl = Ex(l) = a El- Then we have 

i([X,X L ])u L = Cx{i{X l )ujl) - i{X l ){C x ujl) = C x {dE L ) - t(X L )(aio L ) = . 

Hence [X, Xl] is in the kernel of ul and, as ul is symplectic (and therefore regular), the 
kernel is trivial and we arrive at [X, Xl] = 0. 
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(ii) Since X commutes with X L we then have Cx L (i(X)6 L ) = i(X)(Cx L QL)- Therefore 
£ XL (i(X)0 L ) =i(X)(CxM = i(X)dL = X(L)=aL + dh = aL + X L (h), 
and from here we get 

X L {G)=X L (i{X)6 L -h)=aL. 
The result of this theorem can be used to obtain the following virial type relation 



i- [G(t 2 ) - G( tl )} = -Z— f 2 X L (G) dt = f 2 L dt , 



and consequently, in a periodic motion of period T we can take t\ = and t 2 = T an we 
obtain ((L)) = 0. In the more general case, 



1 



-[G(t 2 )-G(t 1 )} = a((L)), 



t 2 — ti 

where L is averaged in t e [ti,t 2 ] and we can also take the limit when t\ — > —00, t 2 — >■ 00. 

As an example we consider the vector field X that is the complete lift of the vector field 
Y generating dilations in Q, given by 

x=Y<=f; 



i=i 



; 9 , (9 

<9x J dv l 



Then the action of X on the Lagrangian L of the harmonic oscillator is given by 

I(L)=l(imv 2 -imw 2 r 2 j = 2L, 

so that X is a symmetry of dynamical vector field of the harmonic oscillator with a = 2. 
Hence we have the following property 



£^2 — ti 



E 



.1 j 



dL_ 



*2 



(8) 



where on the right hand side L is averaged in £ e [£i,£2]- Once again we can particularize 
the time interval t 2 — 1\ for a periodic T of the system (this result is also true for the limit 
£2 going to infinity since the motion is always bounded and periodic). 

This example can be obtained as a particular case when starting with the general one- 
dimensional case of a natural Lagrangian L given by 

L(x, v) = ^ m v 2 — V(x) . 
Let Y be the vector field in K. given by 



(9) 



whose complete lift X is given by 



X(x,v) = Y c (x,v) = ^(x)— + (y 



(10) 



dx ) dv ' 



then we have 



(X(L))(x,v) = —£(x) V'(x) + mv£'(x) v, 



and therefore the condition Y C L — a L, with a e R is written 



r 2e'(x) = a, 

\ V'(x) = aV(x) 



from where we obtain 



£(x) = -ax + C, 



V{x)=C 2 {x + C 1 )\ 



and we recover as a particular case for a = 2 the above mentioned harmonic oscillator. 

4 The virial theorem for position dependent mass sys- 
tems 

The study of systems with position-dependent mass is receiving a lot of attention recently, 
with the non-commutativity of mass and momentum being a difficulty to be taken into 
account in the quantization process [351 GEJ E7] . From the classical point of view the situation 
corresponds to systems for which the kinetic energy is defined from a non-Euclidean metric. 
Therefore the systems are not invariant under dilation and one should look for alternative 
functions giving rise to a virial-like theorem. The more general theory we have developed 
here allows us to deal with this more general case. 

We now consider the one-dimensional case of a Lagrangian L given by 



where m(x) is a positive differentiable function. 

According to the formalism presented in the preceding section, we must look for a vector 
field Y in R given by and its corresponding complete lift X = Y c ffTUl) such that X(L) = 
a L, with a G R. Then as we have 






the mentioned condition leads to 



2£'(z) + £(x) n(x) = a, 
£(x) V'(x) = aV{x) 
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where = m'(x)/m(x). 
The first equation is an inhomogeneous linear differential equation with a general solution 

= -jfUf ( Cl + %f ■ ( 12 ) 

y/m(x) \ ^ Jo J 

Consequently, the vector field Y is a sum of two vector fields. The first one Y\, which is 
obtained putting a = in the preceding expression, preserves the kinetic term T and it 
therefore represents a Killing vector of the associated metric. The second one Yi is given by 



2 y/m{x) \Jo J dx 

The potential energy V(x) is then given by 



V(x)^C 2 exp(J X ^-dcY (13) 



while the function G is 



GO, v) = m(x) £(x)v= yfa(xj [Ci + ^j VHO d(j v. (14) 

and the virial theorem provides the relation 

((X L (G))) = a((L)). 

As an example we can consider the differential equation for a 1-dimensional nonlinear 
oscillator studied in 1974 by Mathews and Lakshmanan |38j 

{l + Xq 2 )q~Xqq 2 + a 2 q = 0, A>0. 

The general solution takes the form q(t) = Asm(ut + 0) , with the following additional 
restriction linking frequency and amplitude 



X 



2 « 2 



1 + A A 2 ' 

The system admits a Lagrangian formulation with Lagrangian: 

L\(q, q) = \ : \ 2 if - a 2 Q 2 ) ■ 
2 1 + A q 2 

It is a system with nonlinear oscillations with a frequency (or period) showing amplitude 
dependence. We can also allow negative values for A, but when A < the values of x are 
limited by the condition \x\ < 1/a/|A| [39]. 



10 



In the limit A — > we recover the equation of motion and the Lagrangian of the harmonic 
oscillator and the frequency becomes independent of the amplitude. In this sense the system 
can be seen as a deformation of the harmonic oscillator. It can also be considered as an 
oscillator with a position-dependent effective mass which depends on A, m\(q) = (1 + A g 2 ) -1 . 
A quantum version of this model was studied in , and a superintegrable generalization to 
several degrees of freedom in classical mechanics was proposed in [10] and the corresponding 
quantum version was studied in |41]-|43j. 

In order to deal simultaneously with all possible values of A it is convenient to introduce 
the following K-trigonometric functions, where k G M: 

cosa/ki if k > 0, ( -^siu^/kx if ac> 0, 

<::„■(.'■) = { i if « = o, s K (x) = I x if k = o, 

cosh-v/— k x if k < 0, y smh-^J—Kx if k < 0, 

and the K-dependent tangent function T K (x) defined in the natural way, T K (x) = S K (x)/C K (x). 
The fundamental properties of these curvature-dependent trigonometric-hyperbolic functions 
are 

C«(x) + kS«(x) = 1, 

and 

C K (2x) = Cl(x) - k S 2 K {x) , S K (2x) = 2 S«(x) C K {x) . 
It can also verified that the derivatives of these functions are given by 

-j-Sk(x) = C K (x) , ^-C K {x) = -k$ k (x) , 

ax ax 

as well as 

d 1 d , 1 



In the case we are considering, the position dependent mass m(q) is given by m(q) 
[1 + Ag 2 ) -1 and the function £(g) becomes 



dq' 



+ Xq' z 

Under the change of variables q = S K {u) with k = —A, and therefore a/1 + A q 2 = C K ( U ), we 
obtain 

e(«) = c(u) (ci + |«). 

or rewritten in terms of the original variables 
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The potential energy V(q) is then given by the integral 



' a ,\ „ / / " a 



that leads to the value 

^)=c( 1+ ^ M y=c(i + ^s- (9 )) 2 . 

For such a potential the function G is given by 

Gfet,) = 7m?( Ci + i s - 1M )' 

and the virial theorem provides the relation 

((X L (G))) = a((L)). 

We close this section with a comment on the general case of the Lagrangian ( JTTT) . If we 
introduce a new coordinate u given by the relation 

u(x) = / \/m(C) d( 



and such that du/dx = ym(x), then the expression for the function becomes 



\Jm{x) V 2 

and the potential V(x) is just a quadratic function in the function u(x) 

V(x)=C 2 (l + ^-u(x)y . 

5 The virial theorem and non-strictly canonical trans- 
formations 

In classical mechanics a transformation that preserves the Poisson brackets up to a multi- 
plicative constant, called the valence of the transformation, is said to be a canonical trans- 
formation [H] (if the valence is the unity the transformation is strictly canonical, otherwise 
is non-strictly canonical). In differential geometric terms, a vector field X on a symplectic 
manifold (M, u) is the generator of a one-parameter group of non-strictly canonical trans- 
formations if there exists a real number a^O such that Cxoj = aw. 
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Given such a vector field X, let us choose a vector field X\ such that 

£ Xl u = -u. (15) 

For instance, when (M, u = —d6) is an exact symplectic manifold the vector field Xi can be 
chosen to be defined by (see [25] ) 

i(X{)u = 9, 

because then 

C Xl 6 = i{X x )de + d{i{X 1 )9) = -i(X x )u + d{i(X x )6) = -9 + d{i{X 1 )9), 
and therefore 

Lx\U = — £-Xi(d9) = —dC Xl 9 = d9 = —ui. 

Note, however, that manifolds exist that do not admit exact symplectic forms, because of 
topological obstructions. For instance, there is no exact volume form in a compact manifold 
M, as a consequence of Stokes theorem, while if u is an exact symplectic form in a 2n- 
dimensional manifold M, then u An would be an exact volume form in M. 

Anyway, if Xi satisfies (IT5|) . then X + aXi is a locally-Hamiltonian vector field since it 
satisfies 

Cx+ax^ = au — au = 0. 
That means the existence of a closed 1-form a such that 

i(X)cu + ai(Xi) u = a. 

Conversely, given a closed 1-form a the preceding relation defines a vector field generating 
a one-parameter (local) group of non-strictly canonical transformations <p t with valence e ae . 

In a Darboux chart for which 

n 

uj = dq l A dpi , 
i=l 

we can choose as 9 (such that u = —d9) the 1-form given by 

1 n 

6 = 2 fa dqt ~ ql dpi ) ' 
i=i 

and then the vector field X\ is the dilation generator 




If a = d(f) (at least locally), then we can write X as follows 

X = Xj, — a, X\ , 
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and then the action of X on the Hamiltonian is given by 

X(H) = {H,<f>}-aX 1 (H). 

Using a Darboux chart, if we assume that H(q,p) = H (p) + V(q), with H a quadratic 
function, we find 



and taking into account that 



n 

a 



we obtain 

X{H) = ^H,<P + ^J2 <?Pi j + 2 a H ° 

Now integrating in time from t = to t = T and, as according to (EJ the average of the 
Poisson bracket vanishes, we obtain, in the limit of T going to infinity (or when the motion 
is periodic with periodic T), the following equality 

2a((H )) = ((X(H))). 

Finally, if there exists b G R such that X(H) = bH, then {H, 0} — aXi(H) = b H, and we 
arrive at 

2a((H )} = bE. 

As a concrete example, let us consider a homogeneous potential of degree d ^ 2, i.e. 



fc=i 



and the one-parameter group of transformations generated by 



n 



1=1 



a — l s d a + 1 9 



Pi 



2 * dqt 2 r dp 



One immediately computes 

£x a w = aw 

and, therefore, X a generates non-strictly canonical transformations. Besides, if we chose 
a = (d + 2) /(d - 2) we have 

Xa (#) = 
v ; d-2 
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i.e. b = Id I id — 2), and we are in the situation described before. Hence, using the previous 
result we have 

(d + 2)((H )) = dE, 

or equivalently 

2((H )) = d((V)), 

which is the standard form of the virial theorem for homogeneous potentials. Another 
example combining scaling of momenta with translations in coordinates was used for the 
Toda lattice in [28J . 

Finally, it is to be remarked that the case we considered in the preceding section of in- 
finitesimal groups of point transformations such that X(L) = aL (up to a gauge term) is 
only a particular case (with a = b) of the preceding situation for the Hamiltonian system 
(TQ,u L ,E L ). 

6 The virial theorem in quantum mechanics 

There has been an increasing interest in the geometric formulation of quantum mechanics 
as a particular case of a Hamiltonian dynamical system (see e.g. [J5|-[I7] and the references 
therein), the main difference with the classical mechanics case being that the symplectic 
manifold is a (may be infinite-dimensional) Hilbert space seen as a real Banach space. More 
specifically, a separable complex Hilbert space ("H, (•, •)) can be considered as a real linear 
space, then denoted "Hr. The norm in "H defines a norm in "Hr, where |H|r = \\v ||c- Of 
course, more accurately one knows that pure states in quantum mechanics are rays rather 
than vectors, but this fact can also be taken into account by replacing the Hilbert space for 
the corresponding projective Hilbert space. 

The real linear space Hr is endowed with a natural symplectic structure as follows: 

cj(u, v) = 2 Im (u, v). 

In fact, a; is a skew- symmetric real bilinear map and the M-linear map Q : Hr — > defined 
by uj{u)v = u(u,v) is not only injective but it is also an isomorphism, because if oj{u) = 0, 
then Q{u){iu) = 2 (u,u) = 0, and consequently u = 0. Riesz theorem can be used to prove 
that the map is also surjective. 

The Hilbert Kr can be considered as a real manifold modelled by a Banach space admitting 
a global chart. Moreover, for each v G Hr the tangent space T V T-Lu is canonically isomorphic 
to the own Hr: we associate with w G Hr the vector in the tangent space T v l-i^ defined by 

Xv{w)f = ^f(v + tw)\ t=0 , 

where the function / is differentiable in a neighbourhood of the vector v, f G C°°(v). This 
is an isomorphism Xv '■ %R — > T v 7i K providing us an identification which corresponds to the 
one given by the free transitive action of the Abelian group of translations on H R . 
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One can prove (see later on) that the constant symplectic structure oj in "Hr, considered 
as a Banach manifold, is exact, i.e., there exists a 1-form 6 G /\ 1 ("Hm) such that uj = —d9. 
Such a 1-form 9 G f\^(H) is, for instance, the one defined by 

0( v )\Xv(w)] = -lm{v,w), 

because then u = —d6 is a symplectic 2-form such that 

u(v)(xv(u),Xv(w)) = 21m (u,iu). 

A continuous vector field in % R is a continuous map X : "% K — >■ "H R . For instance, for each 
v G "H, the constant vector field X v defined by 

X v {w) = Xw(v), 

is the generator of the one-parameter subgroup of transformations of "Hr given by 

$(t,w) —w + tv, 
i.e. with the natural identification of T"H K with "Hr x "H K , 

X„ : w (->■ t>) . 

The values at a point of such vector fields generate the tangent space at the point. 
Similarly, for each vector v G "Hr there is a constant 1-form a v in % R , given by 

a„ : w !->■ (t>, u>) . 

Obviously, 

cv(X„ 2 ) = (ui,u 2 ) , 

and therefore 

a Vl+ \v 2 = a v! + Aa^a , VA G K . 
The 1-form 9 defined above satisfies 

9(X v ) = -lm(-,v), 
because according to the definition of the 1-form 9, 

[9(X v )](w) = 9(w)[X v (w)} = 9(w)[ Xw (v)} = -Im (w,v) . 
One can see that X W [9(X V )] takes a constant value: 

X w [9(X v )](u) = [9(X v )(u + tw)] lt=0 = -Im^-(u + tw,v)\ t=0 = -Im (w,v) . 
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This allows us to check that co = —d6, because for any pair v,w G %, as X v and X w commute, 
[X v , X w \ = 0, we have 

-d9(X v , X w ) = -X v 6{X W ) + X W 9(X V ) = -2 Im (w, v) = u(X v , X w ) . 

As another particular example of vector field, consider the vector field X A defined by the 
C-linear map A : T-L — > %, and in particular when A is self-adjoint. With the natural 
identification of T"Hr pa % R x H R , X A is given by 

X A : v ^ (v, Av) G H x U . 

When A — I the vector field X/ is the Liouville generator of dilations along the fibres, 
usually denoted A = Xi and given by A(f) = (v,v). 

Given a self-adjoint operator A in "H we can define a real function in "Hr by 

a(v) = (v, Av) , 
i.e. the evaluation map, which can be rewritten 

a=(A,X A ). 

Then, 

da v (w) = —a(v + tw)\ t=0 = — [{v + tw, A(v + tw))] lt=0 

= 2 Re {w, Av) = 2 Im {—iAv, w) = u(—iAv, w). 

If we recall that the Hamiltonian vector field defined by the function a is such that for each 
w G T V V. = H, 

da v (w) = ou(X a (v),w) , 

we see that 

X a (v) = -iAv. 

Therefore if A is the Hamiltonian H of a quantum system, the Schrodinger equation describ- 
ing time-evolution plays the role of 'Hamilton equations' for the Hamiltonian dynamical 
system (Ji,u,h), where h(v) = (v,Hv): the integral curves of X h satisfy 

v = Xhiy) = —i Hv . 

The real functions a(y) = (v, Av) and b(v) = (v, Av) corresponding to two selfadjoint 
operators A and B satisfy 

{a,b}(v) = -i(v,{A,B]v), 

because 

{a,b}(v) = [cu(X a ,X b )](v) = u v (X a (v),X b (v)) = 2lm(Av,Bv) , 
and taking into account that 

2 Im (Av, Bv) = -i [(Av, Bv) - (Bv, Av)} = -i [(v, ABv) - (v, BAv)} , 
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we find the above result. 
In particular, on the integral curves of the vector field defined by a Hamiltonian H, 

a(v) = {a, h}(v) = —i (v, [A, H]v) , 

that can be rewritten as 

j t (v,Av) = -i(v,[A,H]v), (16) 

and is usually known as the Ehrenfest theorem. This is the starting point for the virial 
theorem in quantum mechanics. 

Note that in the derivation of (fT6l) we have used that Av belongs to the domain of H. 
This is not a restriction if H is defined in the whole Hilbert space, but in many occasions 
the self-adjoint Hamiltonian has a domain that is only dense in TL and if A does not preserve 
this domain, extra boundary terms should be added to the standard Ehrenfest theorem, see 
refs. [IH]-[5T] where the occurrence of this kind of anomaly in quantum mechanics has been 
studied. In this paper we shall assume that we have no boundary effects or, in other words, 
the anomaly is absent and ffTB"]) holds. Also note that identity f fTB"]) does not require that A 
be selfadjoint. 

Finally, for the sake of completeness, we must comment that if we want to consider the 
projective Hilbert space, one can start with the open submanifold TL— {0} of 7i and the action 
of the two-dimensional Lie group C* = C — {0}, which is an Abelian Lie group isomorphic 
to the direct product R+ ® 17(1) of R+ = {A £ R | A > 0} and U(l) = {e itf> \ cp £ R}. The 
corresponding vector fields are and Xf given by 

X d {v) = ^e~ e v {e=0 = -v , 
Xf(v) = -^e~ l£ v\ €=0 = -iv . 



Let \I> be the map, 



*:H-{0}^l + x5H, *00 = ( Nl, iTm ) • ( 17 ) 



where S7i is the subset of unit vectors of 7i 

SU = {v £ 7i | (v,v) = 1}, (18) 

and denote pr : : R + x STL — > R + and pr 2 : R + x STL —> STL the projections of ^(TL — {0}) = 
R + x STL on each factor. The map pr 2 o : v i-> v/\\v \\ identifies (TL — {0})/lR + with STL, 
(TL — {0})/lR + ~ STL, and we obtain in this way a reduced space which can be seen as 
a submanifold jsn '■ STL —> TL, endowed with the 2-form usn — jsn u ^ which is an exact 
{usu = —dOsui with dsn — 3*su®) ^° u ^ degenerate 2-form. 

As a second step in the reduction process as the symplectic action of the Abelian group 
£7(1) = {e l{p | ip £ R}, on TL, preserves the submanifold STL, the set of orbits is STL/U{1) is 
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the projective Hilbert space RH — {% — {0})/C*. Note that the momentum map J : W — > 
u*(l), can be found to be given by 



(J(v),i\) = X(v,v), 

i.e. J(v) = (v, v) and therefore = SH is endowed with a presymplectic form lush that 

is the pull-back of u. Its kernel is generated by the fundamental vector field corresponding 
to il G u(l). Consequently Marsden and Weinstein reduction process leds to a uniquely 
determined reduced symplectic form in the quotient space, the projective Hilbert space. 
Note however that such a symplectic form is not exact, recall that if dimension of T-L is finite 
its projective Hilbert space is compact. Of course the evaluation map a associated to the 
self-adjoint operator A is not projectable and it must be replaced by the expectation value 
function [48J 

(v, Av) 
{v,v} 

which coincides with the evaluation map at points v G SH. The dynamical vector field 
obtained in this way is projectable and gives us the dynamics in the projective Hilbert 
space. 

If the state v is stationary then the above equation ffl6l) becomes an identity since both 
sides are zero. In a generic state, if we integrate between and T we obtain 

(v(T), Av{T)) - MO), Av(0)) = -i f T (v(t), [A, H]v(t)) dt, 

Jo 

and if (v(t), Av(t)) remains bounded, taking the limit when T goes to infinity of the quotient 
of both sides by T we find 

({(v, [A,H]v))) = 0, 

which is the quantum hypervirial theorem when no anomaly, due to boundary terms, is 
present (5TJ. 

Suppose that the Hamiltonian of a quantum system is 

# = Ip. P + F (X), 

where P = — i V and natural units are used. Let now A be given by 

A = % - (X • P + P • X) . 

We first remark that [X • P, H] = [P ■ X, H] because ofX-P-P-X = 3i 
Taking into account the algebraic relation [AB, C] = A[B,C] + [A, C]B one can see that 

[X.P,ip-P] = ^X J [P J ,P.P] + i^[X i ,P-P] P 1 = l P-P, 

i=i t=i 
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while 

3 

[X ■ P, V(X)} = J2 X * i p i> V W] = -i X • W(X), 

i=i 

and therefore 

(v, [X-P,H}v) = t((v,P-Pv)- (v, (X • W(X))r;)) , 
and we obtain the standard quantum virial theorem 

(((,,p.p,»>-(((,,(x.vy(x)H» = o. 

Note that A = | (X ■ P + P • X) is the generator for the dilation subgroup. Recall that if 
a Lie group G acts on M on the left, then if /i is a G-invariant volume, we can define the so 
called quasi- regular unitary representation in (£ 2 (M),/i) as follows: 

(U(g)iP)(x)=iP(g- 1 x). 

When n is not G-invariant but quasi- invariant, in order to get an unitary representation we 
must correct the right-hand side by the square root of the Radon-Nikodym derivative. 

In the case of dilations in the one- dimensional case, M = R, and if G is the dilation group 
there is no invariant measure. The quasi-regular representation then turns out to be 

[U(X)i;](x) = X'^iX^x). 

The expression so defined is a one-parameter group of transformations with canonical pa- 
rameter a such that X = e a , 

[U(a)ip](x) = e- a/2 ^j(e- a x), 

and its generator comes form 



de~ a l 2 ^(e- a x) 



a=\j 



da 

from where we deduce that the generator of this action is 

1 d_ _ 1 ( d_ d_ 

2 dx 2 \ dx dx 

For M = M 3 , the quasi- regular representation is 

[t/(A)^](x)=A- 3 /V(A- 1 x), 

or in terms of the parameter a, 



da 



3 

= -^^( x ) - x • w, 
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and the generator of the action can be written 



3 1 
A=- + x- V = -(x-V + V-x) 

2 2 



In order to better understand the geometric properties of the virial theorem in quantum 
mechanics, it is interesting to study the quantum version of the system with position depen- 
dent mass introduced in Sec. 4. Its quantization has been studied in detail in (39], where 
it was shown that the appropriate description of the system implies the introduction of a 
non trivial metric g = m(x)dx (g> dx in K. Consequently, the Hilbert space is L 2 (M, djj) with 
d[i = yfm{x) dx and then the norm of a function \1/ is given by 

poo 

\\m\\ 2 = I \V(x)\ 2 y/m(x) dx. (19) 



The Hamiltonian operator of the system reads 

H = H + V(x) = ip 2 + V(x) 

where 

P _ i 9 ; 

y/m(x) dx ' 

is symmetric with respect to the scalar product induced by (TT91 and self-adjoint when ap- 
propriate boundary conditions are chosen. Note also that P is a Killing vector of the metric, 
i.e. Cpg = 0, and it is, up to a factor i the generator of the translation group in this 
non-homogeneous space. 

In order to apply the virial theorem we introduce the generator of the dilation group 
A = £(x)d x with the property 

[P,A] = a -P. 

The general solution was obtained in ( 1T2|) and is given by 

ax) = -jL= (Cx + \ [ X VMt)d() . (20) 
^m(x) V 1 Jo J 

Notice that the name of dilation group is justified because of the commutation relation of 
its generator with that of translations and also because of the action on the metric, i.e. 

£a9 = ag, 

that corresponds to a conformal transformation. 
We consider now the potential 

V(x) = C 2 e X p(- f ^y rf c), (21) 
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that verifies 

[V,A] = bV. 

With all these ingredients we have 

(v, [H, A] v) = a(v, H v) + b(v, V v) , 

and consequently 

a(((v,H v))) + b(((v,Vv))) = 0. 

Finally, some comments on Fock approach to the virial theorem in quantum mechanics 
|52j : Starting with an arbitrary wave function we consider the one-parameter family of 
trial functions {0 A = U(X)<f> | A G where U is a linear representation of the dilation 

group. 

When d — 1 the expectation value of the kinetic term H$ is homogeneous of degree —2 
and the potential V is assumed to be homogeneous of degree k: 

(0 A , H <p x ) = X~ 2 (0, H <i>), (0\ V<P X ) = \ k (0, V<f>), 

therefore 

E x = (0 A , H<p x ) = \- 2 (0, #o0> + A fc (0, K0). 
The best approach to an eigenvalue in the family will be by a value of A such that 

^ = -2 A" 3 (0, #o0> + kX'- 1 (0, V<f>) = 0. 

In particular, if is actually an eigenvector, then the extremal is found for A = 1, 

2(<f>,H 4>)=k(<f>,V4>), 
and we reobtain in this way the virial theorem for eigenstates of H. 

7 Final comments 

The standard theory of virial theorem has been revisited from the perspective of the theory 
of Hamiltonian systems in symplectic manifolds. This allows us to consider the classical 
counterpart of the usually called hypervirial theorems and clarify that the theory is not 
related in the general case with the group of scale transformations. In this way we have 
found room for dealing with systems whose configuration spaces are not linear spaces. The 
geometric formalism is valid not only in the Hamiltonian approach (cotangent bundle with 
the canonical two-form) but also in the Lagrangian case (tangent bundle with a L-dependent 
symplectic form); we have shown the usefulness of some infinitesimal point transformations 
that are not symmetries of the Lagrangian for establishing virial-like relations. The example 
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of position-dependent mass has been used to illustrate the theory and the particular example 
of a nonlinear oscillator is explicitly studied. 

The use of the modern symplectic approach to quantum mechanics enables us to prove in 
the second part of the paper that the virial theorem in quantum mechanics appears in full 
similarity with the analogous classical case. We only need to take into account the point of 
view of symplectic geometry in infinite dimensional linear spaces. We hope this analogy will 
be helpful for dealing with various quantum cases. 
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